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ABSTRACT 

Observations from the last decade have indicated the existence of a general class of superluminous 
supernovae (SLSNc), in which the peak luminosity exceeds 10 44 erg s" 1 . Here we focus on a subclass 
of these events, where the light curve is also tens of days wide, so the total radiated energy is order 
10 51 erg. If the origin of these SLSNc is a core-collapse-driven explosion of a massive star, then the 
mechanism which converts the explosion energy into radiation must be very efficient (much more than 
in typical core collapse SNe, where this efficiency is of order one percent). We examine the scenario 
where the radiated luminosity is due to efficient conversion of kinetic energy of the ejected stellar 
envelope into radiation by interaction with an optically thick, pre-existing circumstcllar material 
(CSM), presumably the product of a steady wind from the progenitor. We base the analysis on a 
simple, numerically solved, hydrodynamic diffusion model, which allows us to identify the qualitative 
behavior of the observable light curves, and to relate them to the parameters of the wind. We 
specifically show that a wide and superluminous supernova requires the mass of the relevant wind 
material to be comparable to that of ejected material from the exploding progenitor. We find the 
wind parameters which explain the peak luminosity and width of the bolometric light curves of three 
particular SLSNe, namely, SN 2005ap, SN 2006gy, and SN 2010gx, and show that they are best fitted 
with a wind that extends to a radius of order 10 15 cm. These results serve as an additional indication 
that at least some SLSNe are powered by interaction of the ejected material with a steady wind of 
similar mass. 

Subject headings: circumstcllar matter — shock waves — supernovae: general — supernovae: individual 
(SN 2005ap, SN 2006gy, SN 2010gx) 
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The discovery of extremely luminous transients in 
the last years led to their c lassification as superlumi- 
nous s upernova e, or SLSNe dQuimbv et al.1 [2007L I20TI 
iPastorello et~aDl2010l ISmith et al.ll2010HGal-Yamll2oT^ . 



In many of these events, the light curve is tens of days 
wide, so the total radiated energy is of order ~ 10 51 
erg. The radius of the photosphere at peak luminosity is 
i? p h > 10 15 cm, as inferred from the observed tempera- 
ture, assuming blackbody emission. These high values of 
radiated energy challenge our understanding of the en- 
ergy source and the conversion mechanism of available 
energy to radiated energy. 

The standard model of core collapse supernovae is that 
the explosion initiates a shock wave which propagates 
through the progenitor and deposits about half of the 
explosion energy as thermal energy and half as kinetic en- 
ergy of the ejecta (for a strong shock in an ideal gas). Red 
supergiant (RSG) supernovae progenitors have a typi- 



cal initial radius of i?„ 



10 



13 



10 14 cm. For such a 



progenitor, as the shock reaches the stellar surface the 
star is highly opaque, so most of the thermal energy is 
transformed to kinetic energy of the ejecta via adiabatic 
expansion, and the fraction of thermal energy which es- 
capes as radiation when the exp anding ejecta becomes 
transparent is roughly i?*/i? p h (|Arnettl Il996f ). Other 
progenitors are even more compact. Therefore, in or- 
der to explain the radiated energy in SLSNe, we must 
either assume an energy source much larger than ~ 10 51 
erg (typical for core collapse SNe), or find a more effi- 
cient mechanism of transforming the explosion energy to 
radiation. 



There have been several suggestions of efficient mecha- 
nisms for conv erting the explosion ene r gy into emerg- 
ing radiation (ISmith fc McCravl l2007t IWoosley et all 



l2007tlChevalier fc IrwiiJ20lHIChatzopoulos et al 
These suggestions are all based upon extensive pre- 
explosion mass-loss - a wind - by the progenitor star 
which expels matter to large radii, comparable with i? p h- 
In principle this can be the result of a burst of mass loss 
or of a steady wind. This matter contributes to the ra- 
diated energy in two complimentary manners. The ex- 
tended mass defines an effective radius R% > i?*, over 
which the shock deposits the explosion energy (essen- 
tially, creating a "bloated" star). With this larger ra- 
dius, less of the thermal energy is lost through adiabatic 
expansion, and a larger fraction i?° ff /i? p h of this energy 
emerges as radiation. In addition, this material interacts 
with the expanding ejecta. This interaction converts ki- 
netic energy of the expanding ejecta back to thermal en- 
ergy via shocks propagating forward through the wind 
mate rial and backward through the ejecta (|Gal-Yaml 
IMl . In the roughest approximation, this interaction 
between the ejecta and the wind can be regarded as a 
plastic collision between an ejecta mass M e j moving at 
a certain velocity and a "stationary" wind mass M w . In 
this approximation, a fraction ~ M W J (M j + M w ) of the 
ejecta kinetic energy will be converted to shock energy. 
If these shocks deposit the energy at a radius ~ R ph 
then it will escape as radiation without suffering signifi- 
cant adiabatic losses. We conclude that an excess mass 
of M w ~ M e j extending to a radius ~ R ph can con- 
vert a substantial portion of the explosion energy to ra- 
diated energy and therefore explain the high luminosities 
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of SLSNc with a conventional ~ 10 51 erg energy source. 

Most of the previ ous works on supernova explosions in 
a den s e mass loss (ISmith fe McCravll2007t ISmith et al l 



2010t iBalberg fc Loebl 120111: iChevalier fe Irwinl l201lF 
Chatzopoulos et al.l 120121 ) used simple approximations 
which allow analytical order of magnitude estimations. 
Some of these wor ks also drew conc lusion from the self- 
similar solutions of Chevalier (1982), which describe the 
interaction of an expanding ejecta with a wind in the 
limit where the wind mass is small, M e j 3> M w , so only 
the outer layers of the ejecta interact with the wind, 
and M e j does not create a natural scale in the dynam- 
ics. For a more accurate relation between parameters 
of the star- wind system and the observed light curve, 
a self-consistent hydrodynamic calculation is required. 
Early numerical calculations of supernov a explosion into 
a den s e win d have been carried out by iFalk fc Arnett 
l973 LfT97l . and more systematically by iMoriva et al 



20111 ) . These works do not cover well the regime M u 



M e j which is of interest in the current work, in order to 
achieve a high efficiency in terms of generating a lumi- 
no us light curve . One exception is the very recent work 
of IMoriva et al.1 (|2012l ). who conducted several numeri- 
cal calculations of a light curve for the specific case of 
SN 2006gy: they considered the relevant mass regime for 
a similar but different case of a shell set at a standoff 
distance i? p h from an exploding star. 

In principle the circumstellar material (CSM) can be 
the result of a steady wind or a burst of mass loss during 
the last stages of the stars evolution. When the wind is 
steady the final configuration is a star enclosed in a con- 
tinuous envelope, while a burst ends with a thin shell of 
matter situated at some distance from the star. In this 
work we focus on the steady wind scenario, for which 
the density profile is well defined (r~ 2 dependence), and 
can be motivated b y various aspects of late stages of 
stellar evolution (see iQuataert fc Shiodell2012t ). We use 
hydrodynamic diffusion calculations to conduct a general 
survey of the relations between the light curve and the 
progenitor and wind parameters. We focus on the par- 
ticular case of a steady wind generated during the last 
stages of the progenitor star's evolution, so the wind es- 
sentially extends from the surface of the star. Our goal is 
to keep the model as simple as possible, and therefore we 
adopt many simplifying approximations, which while still 
sufficient to investigate the key features of light curves 
in the interacting ejecta and wind scenario, are also sim- 
ple enough to be instructive and the results are easily 
understood. The simplifications we applied are detailed 
below. 

The outline of the paper is as follows. In §[5] we present 
our model for general SLSN progenitor systems. In § [3] 
we consider the qualitative relations between the sys- 
tem parameters and the observed light curve, and an- 
alytically examine various limits of the problem. The 
hydro-diffusion code is described in § [4] and numerical 
results are presented in § [5] The observed light curves 
of three specific candidates for an interacting ejecta and 
wind scenario, SN 2005ap, SN 2006gy, and SN 2010gx 
are analyzed and reproduced in § [6j We summarize our 
conclusions in § [7J 



2. PROGENITOR SYSTEMS 



In our scenario the progenitor systems is composed of 
two parts: the gravitationally bound compact star, and 
the outer material which is created by pre-explosion mass 
loss. Both are characterized by their total masses and re- 
spective density profiles. Since the supernova explosion 
induces temperatures much higher than the initial star 
temperatures, and ejecta velocities much higher than the 
initial wind velocity, both the star and wind are approx- 
imated as initially stationary and cold. The explosion 
can then be treated as an instantaneous release of ther- 
mal energy E at the center of the star (r = 0). 

For systems in which the photosphere lies at R p h ^> R* 
and M j ~ M w the details of the structure of the pro- 
genitor star are unimportant in terms of the resulting 
light curve, since the star can be considered as a point 
object. None the less, we do want to allow also for 
calculations of systems where the wind mass is much 
smaller than the ejected mass, mainly for comparison 
with other works w hich invest igated these cases (see e 
IMoriva et aI1l2011l) . We follow lMatzner fc McKeel (| l99 ; 
by modeling the progenitor star by a polytrope with ra- 
dius i?*, mass M*, and polytropic index n. n = 3/2 
is suitable for convective envelopes, like those in RSGs, 
while n = 3 is used to describe radiative envelopes, such 
as those in blue supergiants (BSGs). We approximate the 
star by applying the Lane-Emden equation to the whole 
star, not only to the outer envelope. Thus, i?*,M*, and 
n completely determine the initial density profile of the 
star. Note that since the Lane-Emden equation is ap- 
plicable only for the envelope, M» actually denotes the 
envelope mass, and not the star mass. All these ap- 
proximations are valid since the inner structure of the 
mantle has little effect on the shock wave propagation. 
The mantle can be co nsidered as a point mass which is 
added artificially (see iMatzner fc McKeel Il999l) . Thus, 
in terms of size, a 15Mq progenitor should actually be 
modeled by a M* ~ IOMq structure, to allow for the 
~ 5Mq mantle. 

For the wi nd density profile, we followlBalberg fc Loebl 
(|2011l) and IChevalier fc Irwinl (|2011[ ) and consider a 

steady wind with mass loss rate M and wind velocity 
v w . The resulting density profile is 



p{r) 



M 



Kr~ 



(1) 



We impose this wind density profile by continuously at- 
taching it to the polytropic structure of the star at a 
radius where the polytrope and equation ([1]) coincide. 
For the parameters chosen in this work, this wind base 
radius is essentially at 

The wind is assumed to continue to an outer radius 
R w so that its total mass is 



Ai:r 2 Kr 2 dr 



AttK{R w -R,). 



(2) 



For R w 3> i?* , which is the case of interest in this work, 
equation ([2]) is reduced to 



M w w 4nKR u 



(3) 



For numerical convenience, we smooth the wind edge at 
R w . Instead of cutting the wind abruptly at R w , we 
extend the wind to a larger radius r — uR w with a > 1 
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(we usually choose a ~ 1.2) and change the wind density 
profile from equation (fT]) to 



p(r) = Kr 2 e 



(4) 



with k ~ 10 providing a sharp power law. The details of 
this smoothing have a minor effect on the results. 

Since the light curve is mostly dominated by the struc- 
ture of the wind material, we mostly survey the wind pa- 
rameters and fix the parameters of the star. We choose 
a star with R* = 10 13 cm, M„ = 15M Q , and n = 3/2. 
These values are representative of a typical RSG. We 
study the case R w » i?*, for which the value of i?* has 
little effect on the light curve. The value of n (3/2 or 3) 
also does not effect the results in the regime of interest. 
However, the value of M* does has a significant effect 
on the light curve, and the choice of 15M is arbitrary. 
Some consequences of changing M* are discussed in § (6J 

The motivations described in § [1] lead to the choice 
R w ~ i? p h <~ 10 15 cm (the relation R w ~ i? p h is justified 
in §[3]) and M w ~ M e j ~ M*, which, using equation ([3]), 
sets K ~ 10 18 g cm -1 . 

To summarize, our progenitor system is thus modeled 
by 4 parameters: the star mass M*, the explosion en- 
ergy E, the wind outer radius R w , and the wind density 
coefficient K from equation (TTJ). After fixing the star 
mass, we are left with 3 parameters which determine the 
properties of the light curve. 

3. QUALITATIVE PICTURE AND EXTREME LIMITS 

The supernova explosion initiates by a strong shock 
wave, which when close to the edge of the star becomes 
radiation dominated. Diffusion of energy carried by ra- 
diation causes the shock to devel op a finite wid th with 
optical depth St ~ = c/u s h (| Weaver! 1 1976|) . where 
u s h is the shock velocity and c the speed of light. This op- 
tical depth can be intuitively understood by equating the 
hydrodynamical time scale and the diffusion time scale 
over the shock front. At large optical depths from the 
surface, r 3> /3T , the shock wave can still be treated 
as an ideal discontinuity, and diffusion can be neglected. 
When the shock approaches the surface, and r ~ 
energy can escape by diffusion to the surface, and and 
the shock dissolves; the purely hydrodynamical (with dif- 
fusion neglected) solutions are no longer valid. At this 
stage the shock is said to "breakout". 

Shock breakout through the outer layers of a bare star 
has been discussed in severa l works. These rely on the 
self similar solution found bv lSakural (|1960f ) for a planar 
shock propagating through the steep ly declining density 
profile at the edge of the star; see Matzncr fe McKed 
(|1999l ) for a complete derivation. ISapir et al.l (|2011| ) ex- 
panded this self similar solution to include diffusion and 
thus calculate a light curve under these conditions, again 
- for the plane-parallel case. In the presence of an ex- 
tended wind, due to a pre-cxplosion mass loss from the 
progenitor, shock breakout must be considered in spher- 
ical symmetry. As the outgoing ejecta plows through the 
wind and slows down, it drives a forward shock through 
the wind and a reverse shock that propagates back into 
the ejecta, and it is the combined shock profile which 
eventually breaks out as it reaches a low optical depth 
region in the wind. Our focus is on the case where 
M w ~ M c j, so the entire ejecta participates in the shock 



breakout, and the light cur ve is the result of the bre akout 
through the opaque wind (jChevalier fe Irwiiil 120 111 ). 

While no self-similar solution exists for the general 
case, we can still quantify several important conclusions 
about the typical time scale of the light curve. We can 
also draw some insight from the two extreme limits where 
M w ;§> M e j and M w <C M e j. If the wind mass is much 
larger than the ejecta mass, the dynamics are similar to 
the Sedov- Taylor explosion, but with a power law am- 
bient density instead of a constant one. In the opposite 
extreme l i mit (M w <C Af j), the self-similar solution of 
IChevalierl ((19821) can be used to model the combined for- 
ward and reverse shock. For such a configuration only 
the initial light curv e will be the result of s hock breakout 
through the wind (jBalberg fc Loebl 1201 ll ). followed by 
the main light curve driven by the internal energy held 
by the bulk of the ejecta. In the following we examine 
the estimates that can be made about the resulting light 
curve from these solutions and other considerations; the 
qualitative picture we draw serves to clarify the numeri- 
cal results presented in the later sections. 

3.1. Breakout Radius 

As the forward shock propagates through the wind, 
it breaks out when it reaches optical depth r ~ 
c/v s h- We adopt the assumption of constant opac- 
ity, which is appropriate for electron (Thompso n) scat- 
terin g (a similar assumption was adopted by Arnettl 
T996l IChevalier fc Irwint |20TH iBalberg fc Loebl 120111: 
Moriva et ail 120 111) . The value of k is, of course, com- 
position dependent, ranging between k ~ 0.2 cm 2 g _1 for 
Hydrogen free matter and k f» 0.4 cm 2 g _1 for pure Hy- 
drogen (we usually used k w 0.34 cm 2 g _1 which is ap- 
propriate for a 70% Hydrogen composition). In this case 
of constant opacity the optical depth from a radius r to 
the edge of the wind is given by 



r(r) = 



updr = 



nKr 



- 2 dr = nK {-- — 

R w 



(5) 

Therefore, the shock breaks out at radius i? s h which sat- 
isfies 

V s h V ^sh tUu J 



Following Chevali er fc Irwiiil (|2011l ) we denote 

_ nKv sh 
= , 



(7) 



and rewrite equation ([5]) as 

111 

Rsh Rw Rd 

IChevalier fc Irwir] (|2011|) discussed the breakout at the 
two limits: 

R w Rw *C Rd 
Rd R w 3> Rd 



(8) 



(9) 



Our focus is on the intermediate case of Rd ~ R w , which 
appears to be motivated by observations. Note that as 
long as the wind mass is of order the ejecta mass, the 
shock velocity naturally tends to t> s h ~ y/E/M* ~ 3 x 
10 s cm s _1 , so for the parameters chosen in § [21 Rd ~ 
10 15 cm. 
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3.2. Luminosity and Timescale 

The main features of the observed light curve are to- 
tal radiated energy, the peak luminosity and the typical 
width. The three are connected, of course, through the 
emission of the thermal energy at shock breakout by pho- 
ton diffusion through the wind. Being an integral quan- 
tity, the total emitted energy, -E ra d is expected to follow 
the plastic collision picture described in § [1] giving the 
simple relation 

E iad ocE ^ Mw ^ ~E „ M \ r . (10) 



M* + M u 



The typical timescale in the light curve must depend 
not only on integral quantities, but on the details of the 
wind profile as well. We estimate the time it takes the 
shock energy to diffuse to the surface following breakout 
as follows. First, we note that the photons do not diffuse 
all the way from i? s h to the surface, but rather to the 
radius of the photosphere i? p h which is located at r ~ 
1 (specifically r = 2/3 for Eddington's approximation). 
Using equation ([3]), i? p h is given by 

1 ~ 1 + 1 (11) 
-Rph Rw kK 

For the case R w <C Rd < kK, equation (jTTj) yields R p h ~ 
R w . According to equation ([8]) i? s h ~ R w in this case as 
well, so the shock breaks out very close to the edge of the 
wind and we can assume a constant diffusion coefficient 

c cR 2 

D~-~ C -^ (12) 
up kK 

and a diffusion distance (see equations (jHJl, (|11|) . and as- 
suming v s h <C c) of 

B „ „ - R 2 



AR 



Rsh ~ Rw — Rsh 



Ru 



R d 



r>2 

R d 



(13) 



which results in a diffusion time 

, Ai? 2 Rl 
td 



D 



(14) 



In the opposite limit, R w 3> Rd, and so i? s h w Rd <C R w - 
Since the shock velocity is far from rclativistic, we also 
have R s h ~ kKv^/c <C kK. These relations imply, 
according to equation (|TT|) . that i? s h <C R p h- In this 
case shock breakout evolves through the diffusion of ra- 
diation to the photosphere, and the diffusion time can 
be estimated by considering the change of the density- 
dependent diffusion coefficient 

td ps [** d( - T Z Rsh ^ „ [ Rph t ~ R sh)Kpdr , 

jR sh c 
Kr~ 2 density profile the inte- 



Jn* D{r) 

After substituting a p 
gration yields 

kK 
t d « 



hi 



R 



P h 



- 1 



(16) 



Rsh. 

Rsh Rph 

Combining this result with equations flSJ) and (fTTj) gives a 
general expression for the diffusion time for any R w . By 
substituting i? s h ~ Rd from equation (JT)) and i? p h from 
equation (|11|) . and using i? s h <C i? p h we summarize the 
diffusion time at both limits: 
Rl 



hi 



v ah 1 + [kK/R w ] 



1 



R w <C Rd 
R w 3> Rd 



(17) 



Equation (jXTJ) gives a diffusion time which is monotoni- 
cally increasing with R w in both limits, but more mildly 
at large R w . For R w ^> kK ^> Rd the diffusion time 
reaches an asymptotic value of 



kK 

t d ^ 

c 



k/C , / c , 
In — . 18 



The expression i ~ kK/c, sometimes with the 
logarithmic c orrection mentioned, was used in pre- 
vious works (lOfek et all l2010t iBalberg fc Loebl 120111: 
IChevalier fc Irwinl l201l( T to estimate the order of mag- 
nitude of the light curve time scale. Indeed, for K ~ 
10 18 g cm -1 (see § [J), this time scale is ~ 100 days, 
which is cons istent with the time scale in the relevant 
observations (iQuimbv et aT]l2007l l2011t [Pastorello et al.l 
[2010t ISmith et al.ll2010t ). However, since we are inter- 
ested in the regime where Rd and R w are of the same 
order of magnitude, a more careful analysis is necessary 
and we must take into account the dependence of td on 
Rw, which is evident in equation (TT71) . Note that the pa- 
rameters chosen in §[2] do dictate that R w kK, which 
means, according to equation (jlip . that the radius of the 
photosphere i? p h is roughly the wind radius R w . Thus, 
our choice of R w ~ 10 15 cm, motivated by the observed 
photosphere radius (see § [lj is self consistent. 

In general, the shock velocity changes with radius, and 
we can roughly estimate that in the regime M u 



Mb 



Vsh OC 



E 



Mr 



1/2 



E 



M* + AnKR sh 



1/2 



(19) 



with M s h oc KR s h is the accumulated mass enclosed by 
i? s h- When w s h is a function of R s h, and therefore a 
function of R w the functional dependence of td on R w 
can be more complex than presented in equation (|17[) . 
Nonetheless, we can gain significant insight by relateing 
v s h to the progenitor system parameters in the limits 
where self-similar solutions exist. The fundamental point 
is that for a fixed value of K, the limit R w ^> Rd also 
corresponds to M w ^> M c} , while in the opposite case of 
R w <C Rd we also have a low mass wind, M w <C M e j. 

In a very massive wind dim ensional analysis similar to 
the Sedov- Taylor problem fsee lZel'dovich fc Raizer|[l966|) 
can be applied, leading to a shock radius and shock ve- 
locity which evolve as 



-Rsh oc 



Ef_ 
~K 



2 \ 1/3 



Vsh = Rsh OC 



E 



KR 



sh 



1/2 



(20) 



However, since in this limit breakout occurs at i? s h ~ Rd, 
independent of R w , the shock velocity at breakout is also 
independent of R w . Equation (fTTj) provides a qualita- 
tive understanding of the dependence of td on K and 
E; neglecting the logarithmic correction, we have (for 
R w 3> Rd) 

t d oc K l E°. (21) 

For a very low mass wind, v s h can be found wit h the 
alternative self-similar solution of I Chevalier! (|1982h who 
studied the interaction of ejecta, with the unshocked den- 
sity profile p e j(r,t) oc r _ "T" 3 , and wind with density 
profile Pw (r) oc r~ s . In our ca se s = 2 and m can be re- 
lated (see iMatzner fc McKeej[l99l iRabinak fc Waxmanl 
|2011[) to the poly tropic index n: m ~ 10 (for n = 3) or 
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t d OC ^-1^2/(^-3)-! 



m «i 12 (for n = 3/2). The shock propagates in time 
with radius 

i?sh K E^ ( - m ~ 3 ^ ( - m ~ s ' > K~ 1 ^ m ~ 3 h < - Tn ~ 3 ^ < ' m ~ s \ (22) 
and with velocity given by 

V ah = R sh OC Bl/2^-l/C™-3) i ^-3)/(m-3) - (23) 

Combining equation (|23p with equation (|1 T[) results in 
an inverted dependence of the diffusion time on K in the 
regime R w <C i?^: 

'E^K- 5 ' 7 m w 10 (n = 3) 
E -i K -7/9 mPS 12(n = 3/2) 

(24) 

In this regime i? s h ~ -Rw and, using equation (|23[) . we 
have 

Wsh oc fl-V(-3) „ m«10(n = 3) 

\i?^ 1/9 m«12(n = 3/2) V ; 

This is a weak dependence on i?,„, which would result 
in a small deviation from the relation td oc i?^ of equa- 
tion (|T7| . The inversion of the dependence of on K 
between limits (linear in K for R w ^> i?^ and inverse 
for R w <C suggests a weak dependence on K for 
the intermediate regime. We note that, in theory, if the 
Sedov- Taylor solution were applied to the R w <C Rd case, 
by combining equation (|20|) with equation (|17j) , we would 
also have the ^ ~ dependence, implying the robust- 
ness of this result. 

We conclude that when identifying the limits of R w 3> 
R d and R w < R d with Af lu > M ej and M m <C M c j, 
respectively, then at both limits Vsh is approximately 
constant (a different constant for each limit), justifying 
equation (|17|) . We expect the region M w ~ M e j to ex- 
hibit some deviation from these time scale estimates. 



4. THE CODE 

We have written a ID Lagrangian computer program in 
order to calculate the shock propagation and light curve. 
In this section we describe the code in brief. The code 
uses the standard vo n Neumann and Richtmyer stag- 
gered mesh method (Jvon Neumann fe Richtmverl 119501 : 
iRichtmver fe Morton! Il967| ) to solve the nonrelativistic 
equations of motion. The energy equation is solved im- 
plicitly, and the radiative flux is added to the hydrody- 
namics assuming local thermal eq uilibrium (LTE) and 
in th e diffusion approximation (see IZel'dovich fe Raizerl 
|1966[) . More specifically, the code solves the following en- 
ergy equation (implicitly, solving a tridiagonal equation 
system) : 



de 
F 



9V IVF.O 
P 

L»V (oT A ) 



(26a) 
(26b) 

where e is the specific energy, V = 1/p the specific vol- 
ume and p is the pressure. F denotes the radiative flux, 
with D = Ac/ 3 as the diffusion coefficient, and a the ra- 
diation constant. The temperature T in equation (|26b[) 
is the temperature of the fluid (a result of the LTE as- 
sumption). In the case of constant opacity, the mean free 
path satisfies A = 1/ up. Gravitation can be neglected for 
the description of the shock propagation and breakout, 
since GM%/R* < E. 



The equation of state (EOS) is that of a perfect gas, 
with radiation terms added to the pressure and energy: 

p( p ,T)=TpT+^ (27) 



rr 



3 

aT 4 



e(p,T) = —— + — 



(28) 



7-1 p 

where r = R/ p is the gas constant divided by the molar 
mass. In the examples shown in this section we choose 
7 = 5/3, suitable for monoatomic gas and p = 0.6 which 
corresponds to a fully ionized mixture of hydrogen and 
helium with primordial ratios. In the context of shock 
breakout the energy and pressure of the fluid are dom- 
inated by radiation after the shock passage and during 
the period of adiabatic expansion which follows, so these 
terms have a minimal effect on the equation of state and 
the light curve. 

The assumption of thermal equilibrium and the va- 
lidity of the diffusion approximation ar e reasonable for 
nonrelativisitc shocks (jKatz et al.ll2010t ). and we rely on 
them here for our estimates. We do note that near the 
photosphere during shock breakout ( Nakar fc Sarill2010t 
iRabinak fc Waxmanll201ll ) transport, rather than diffu- 
sion, is a more appropriate description of the photon 
propagation, which we will apply in future work. 

As a code check, we calculated several test problems 
with our program and compared the results to known 
solutions. Among the test problems we cons idered are 
the s elf similar interaction of ejecta and wind (|Chevalierl 
Il982h . Elliot's extension to th e Sedoy- Tay lor explosion 
which includes r adiative flux dElliotl Il960h . and planar 
shock breakout (|Sapir et"aT1l2011|) . We present our re- 
sults for the planar shock breakout in the Appendix. 

We present an example of a calculation in figures [T] 
and [3J Figure [T] shows the temperature profiles at differ- 
ent times. The transition from a discrete shock front to 
breakout at t ~ c/w s h, as discussed in § [3j is evident in 
the late time (50 days after explosion) profile. In figure 
[21 the corresponding light curve is presented. It can be 
seen that the emergence of a light curve and its duration 
correspond to the shock breakout time scales. 



4.1. Comparison with \Moriva et all 1(201 A) 

As a fi nal test, we compare d our numerical results with 
those of iMoriva et al.l (|2011[ ) who also studied the effect 
of a wind on light curves. Their work focused on less 
dense winds and using a more complex numerical model 
which we do not try to reproduce. However, we can com- 
pare our calculations to the progenitor-wind model which 
is the most relevant to the regime studied in the current 
work: model sl5w2r20mle3, ma rked as 10 _1 MQyr _1 in 
figures 3,4 of lMoriva et al.l ([201 If ), which has a wind mass 
of 6.5M0. In figure[3]we show a comparison between the 
light curve calculated with our code and their results. 
Our light curve was calc ulated using the prog enitor pro- 
file given in figure 3 of IMoriva et al.l ([20111 ) , which is 
slightly different from our standard A/* = 15Mq progen- 
itor (see § O . The resulting relevant light curve parame- 
ters (peak luminosity, time duration, and total radiated 
energy) are the same up to few percents. 

5. NUMERICAL PARAMETER SURVEY 

In this section we study the light curves calculated with 
the code and conduct a parameter survey in order to 
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Fig. 1. — Calculated temperature profiles at different times fol- 
lowing explosion. The profiles were obtained with the parameters 
R w = 2.5 X 10 15 cm, K = 10 18 g cm" 1 , and E = 5 X 10 51 erg. 
The plotted profiles are from 5 (solid black line), 25 (dashed blue 
line), and 50 (dot-dashed red line) days following the explosion. 
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Fig. 2. — Calculated light curve. The light curve was calculated 
with the same model as figure [l] 



map the dependence on the important parameters. We 
also verify that the analytical limits presented in § [3] are 
recovered, and use them to gain insight to the results. 

5.1. Comparison with analytical limits 

Following a convergence test, we model the compact 
star part of the progenitor system with n* = 50 cells, ge- 
ometrically decreasing in size Ar toward the outer edge 
with a quotient q* = 0.98. The wind is divided into 
n w = 250 cells, geometrically increasing in size (with a 
constant quotient q w > 1) toward the outer wind edge. 
The innermost wind cell is of the same size as the out- 
ermost star cell. The size of the innermost wind cell, 
together with n w , determine q w . In these calculations 
we use the smoothing parameters a — 1.25 and k = 15 
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Fig. 3. — Comparison of the light curv e calculated by our c ode 
(solid black line) to the one calculated by [Moriya et al. ( 20Tlf) for 
the model sl5w2r20mle3 (dashed blue line). For the com parison 
we used the same progenitor profile as Moriya e t alj 11201 If ), 

(see §H]). Generally, we set i?* = 10 13 cm, except for very 
low values of R w , for which we had to reduce the star ra- 
dius, i?*, in order to remain in the regime i?» <C R w - 
Additionally, we had to change the resolution for conver- 
gence and to disable the wind edge smoothing (see § [5} . 
These changes have an effect only at very low R w . 

We relate the diffusion time scale td with the full width 
at half maximum (FWHM) of the light curve. This 
choice is independent of the low luminosity " tail" at late 
times, which arises in part from the continued interac- 
tion of the shock wave and the wind (jChevalier fc Irwml 
1201 lh . and is sensitive to the shape of the wind cutoff pro- 
file. The results, for E = 5 x 10 51 erg, K = 10 18 g cm" 1 , 
and a wide range of R w are shown in figures [4] and [5] 

Both figures demonstrate that the calculated light 
curves do indeed reproduce the analytic timescales at 
the appropriate limits. The fit is very good at low values 
of R w , while at high values there is some deviation from 
equation (|17j) . The reason is that the simplified treat- 
ment of diffusion in the code extends all the way to R w , 
whereas our analytic estimates were based on emission 
from a r = 1 (or r = 2/3) surface. A more exact fit of 
the numerical results is therefore found with a revised 
analytical estimate of the diffusion time using equation 
Ijl6|) . and substituting i? pri = R w : 
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Unlike equation ([T7]). equation (f2T)|) does not reach an 
asymptotic value at R w 3> kK. Correspondingly, our 
numerical calculations which do not take into account 
the photosphere at r <~ 1 are inexact in this sense. How- 
ever, since we deal with the regime R w <C kK, the differ- 
ence between equations (fTT)) and (|29[) is small (sec figure 
[5l which even approaches R w ~ kK), so the numerical 
treatment of the photosphere is sufficient for our pur- 
poses in the current work. 

5.2. Progenitor system - light curve relation 
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Fig. 4.— Full Width at Half Maximum (FWHM) of the cal- 
culated light curve as a function of the wind outer radius R w . 
Results arc^plotted for nominal values of K = 10 18 g cm -1 and 
E = 5 X 10 erg (black line, marked with circles). Each marker 
represents a single hydrodynamic diffusion calcula tion . The an- 
alytical limit (solid blue line), following equation \17). is of the 
form t<i = (kK/c)^ 1 v~ h , with fitted parameter v s ^ = 3.5 X 
10 s cm s _1 
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Fig. 5.— Full Width at Half Maximum (FWHM) of the cal- 
culated light curve as a function of the wind outer radius R w . 
Results are plotted for nominal values of K = 10 18 g cm -1 and 
E = 5 X 10 erg (black line, marked with circles). Each marker 
represents a single hydrodynamic diffusion calculation. The cal- 
culations are the same as in figure [3] but plotted in a different 
scal e. The analytical limit (solid blue line), following equation 
H29t for the approximation R p ^ = R w (which is used by the code; 
see text), is of the form td = r)(nK / c)\\n(R w / R^) — 1], with fit- 
ted parameters 77 = 0.9 and R^ = 1.2 X 10 15 cm (or cquivalently 
^sh = 1.1 x 10 8 cm s _1 ). A second, more adequate, analytical limit 
(dashed red line), with the sam e v s ^ and multiplication factor r), 
is in the form of equation I U7H . This limit takes into account the 
photosphere at r = 1. 



As is obvious from the qualitative discussion in §[3] and 
from figures [4] and [5l the range of R w of interest in this 



work (R w ~ 10 15 cm) cannot fit well to any of the an- 
alytical limits. Henceforth we conduct a numerical pa- 
rameter survey in this range. We begin by examining 
the relation between the diffusion time, td, and the pa- 
rameters of the progenitor system, namely E,K and R Wl 
whereas for simplicity we fix ill* = 15Mq. We study the 
dependence of td on R w for a nominal calculation with 
E = E = 5 x 10 51 erg and K = K = 10 18 g cm -1 . 
The dependence on E and K is studied by repeating 
the calculations with other values of E = E/Eq and 
K = K/Kq. The results are shown in figures [6] and [7] 
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Fig. 6.— Full Width at Half Maximum (FWHM) of the cal- 
culated light curve as a function of the wind outer radius R w . 
Results arc plotted for nominal values K = K/ (10 18 g cm -1 ) = 
E = E/ (5 x 10 51 erg) = 1 (blue line, marked with squares), for 

K = 1/2 (black line, marked with circles), and for K = 2 (red line, 
marked with triangles). Each marker represents a single hydrody- 
namic diffusion calculation. 

The increase of td with increasing R w is understood 
qualitatively by equation (|17|) . Quantitatively, we see 
that td is strongly dependent on R w througho ut the rele- 
vant range, so the approximation td ~ kK/c (|Qfek et all 
[2OT0tlBalberg fc Loebll201ll ; IChevalier fc Irwidl20liF can 
serve only as an order of magnitu de estimate. This be - 
havior is also evident in figure 9 of iMoriva et~ai] (|2011[) . 
As can be seen in Figure [6l the dependence of the diffu- 
sion time on K is in good agreement with the analysis of 
§ 13.21 At large R w , the diffusion time increases roughly 
linearly with K, but for lower R w the dependence be- 
comes weaker and is finally inverted, as expected in low 
mass wind (see equations (|2Tj) . ([M]) ). yielding smaller dif- 
fusion times for larger values of K . We also recover the 
td ~ E~ x relation expected at low R w (again, see equa- 
tion (f24|) ). The dependence becomes weaker at larger 
R Wl as might be expected from equation (|2ip . 

The other feature in the light curve which we relate 
to the parameters of the progenitor system is the total 
radiated energy Figure [5] shows the dependence of the 
radiated energy on K for different values of E and for 
a fixed R w = 2.5 x 10 15 cm. We calculate the radiated 
energy by integrating the light curve until the luminosity 
drops to 0.1% of the maximum luminosity. The behavior 
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Fig. 7.— Full Width at Half Maximum (FWHM) of the cal- 
culated light curve as a function of the wind outer radius R w . 
Results arc plotted for nominal values K = K/ (10 18 g cm -1 ) = 
E = E/ (5 X 10 51 erg) = 1 (blue line, marked with squares), for 

E = 1/2 (black line, marked with circles), and for E = 2 (red line, 
marked with triangles). Each marker represents a single hydrody- 
namic diffusion calculation. 

exhibited in figure [8] is qualitatively understood by the 
plastic collision relation (equation (fT0|) ) . thereby saturat- 
ing when M w > M*. Using ©, M w w A/* for K = 1. 

We digress and discuss the shape of the light curve 
because the total radiated energy is not always reliably 
observed. When observations are limited to the vicinity 
of peak magnitude and do not track low luminosities, the 
shape of the light curve must be modeled theoretically 
to asses its total energy. One such model uses the ob- 
served peak luminosity L max and the observed FWHM 
of the light curve, and, assuming a Gaussian light curve 

(|Arnettj[l99l . £ rad w (o.50r/ln2) L max x FWHM. 

We plot this estimate of the radiated energy (using the 
calculated peak luminosity and the FWHM of the calcu- 
lated light curve) as dashed lines in figured] It is obvious 
from figure [8] that as the wind becomes more dense, this 
simple estimate deviates more from the total radiated 
energy (compare the solid and dashed lines in figure [8|. 
The reason is that the late time "tail" strongly deviates 
from a Gaussian form, and includes much more energy 
than predicted by a Gaussian approximation. This tail is 
partially powere d by continued interactio n of the shock 
wave and wind (j Chevalier fc Irwml |2011| ) and thus be- 
comes more prominent and contains more energy. We 
can conclude that total energy estimates based on obser- 
vations should be done carefully, and that a Gaussian fit 
must be treated as a lower limit if only the vicinity of 
the peak region of the light curve is observed. 

5.3. Time scale constraints 

In this work we focus on winds with radius R w ~ 10 15 
cm. The justification we gave for this choice is the ob- 
served photosphere radius i? p h ~ 10 15 cm and the rela- 
tion R w ~ i? p h which is valid for M w ~ M e j (see § [3]). 
One of the conclusions of the numerical results is that the 
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Fig. 8. — Radiated energy as a function of the wind density coef- 
ficient K, for winds with outer radius R w = 2.5 X 10 15 cm. Results 
are plotted for a nominal value E = E/ (5 X 10 51 erg) = 1 (blue 
lines, marked with squares), for E = 1/2 (black lines, marked with 
circles), and for E = 2 (red lines, marked with triangles). Each 
marker represents a single hydrodynamic diffusion calculation. For 
each value of E the total radiated energy is plotted as a solid line 
and a Gaussian estimate for the radiated energy, based on the peak 
luminosity and FWHM of the light curve (see text), is plotted as 
a dashed line. 



observed light curve time scale gives another constrain on 
R w . We demonstrate this point in Figure where we 
show the light curves calculated for three combinations 
of K and R w , all of which satisfy KR W = 2.5 x 10 33 g, so 
that the total wind mass is kept constant. The intermedi- 
ate model with R w = 2.5 x 10 15 cm is used below to fit SN 
2010gx (see § 0, while the other two are more compact 
and extended winds (note that for the extended wind, 
the approximation R w <C kK which our code assumes is 
marginal). By keeping the total wind mass constant the 
efficiency of converting a given explosion energy to radi- 
ated energy is fixed (for a given progenitor star mass) and 
so the light curve time scale becomes a direct indicator 
of the wind radius. It is obvious that for time scales of 
~ 50 days, as is the case of SN 2010gx, the models with 
the more compact and more extended winds are ruled 
out. 

6. OBSERVATIONS 

In this section we relate our model to SLSNc ob- 
servations. Specifically, we focus on three events: SN 
2010gx, SN 2006gy and SN 2005ap. Our goal is to 
ascertain that a steady wind model can provide a vi- 
able explanation for the observed SLSNe light curves, 
and to constrain the likely parameters of the progeni- 
tor system. In gene ral, previous works which considered 
a sta r -wind system ( Smith fc McCravll2007l : iSmith et al.l 
l2010t I Chevalier fc Irwir] 1201 1[ ) provided only order of 
magnitude correlation between the steady wind mode l 
and the observations. Very recently. iMoriva et al.l ([2012D 
have presented a specific numerical model for SN 2006gy 
where the ejecta interacts with a a distant shell (rather 
than a steady wind) , and we comment on the similarities 
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Fig. 10. — Calculated and observed light curves of SN 2010gx. 
The calculated light curve (solid black line) was calculated with 
the parameters R w = 2.5 X 10 15 cm, K = 10 18 g cm -1 , and E = 
2.5 X 10 33 g constant). The plotted light 2 X 10 51 erg. The obse rved light curve (blue circles) is taken from 



Fig. 9. — Calculated light curves for different wind radii. The 
plotted light curves are obtained from models with the same re- 
leased energy E = 2 X 10 51 erg and the same wind mass (obtained 
by keeping KR U 



curves have wind radii of R w = 5.0 X 10 14 cm (solid black line), 
R w = 2.5 X 10 15 cm (dashed blue line), and R w = 1.0 X 10 16 
cm (dot-dashed red line), corresponding to density coefficients 
values of K = 5 X 10 18 g cm -1 , if = 1.0 X 10 18 g cm -1 , and 
K = 2.5 X 10 17 g cm' 1 respectively. 



IPastorello et al.l l20loh . with zero bolometric correction. 



that our model is compatible with observations.. 



and differences regarding this specific object below. 

6.1. SN 2010gx 

We take the data for SN 2010gx from IPastorello et al.l 
(|2010T) . The measured luminosity and temperature imply 
a blackbody radius ~ 3 x 10 15 cm at peak luminosity. 
Since the photosphere lies close to the wind edge (see 
§ [3]), we choose a model with R w = 2.5 x 10 15 cm (the 
photosphere expands following the passage of the shock). 
We find that the light curve is recovered well when setting 
the other parameters to be K = 10 18 g cm -1 and E = 
2 x 10 51 erg. The quality of the fit is shown in figure ITOl 
which compares the calculated and observed light curves. 
It is noteworthy that we do include a correction for the 
difference in arrival times of photons originating from 
different positions on the photosphere (jKatz et al.ll2012[ ). 
This effect on the time scales is ~ i? p h/c, which is about 
a day, and therefore negligible in the cases considered 
here. 

The corresponding blackbody temperature for the cal- 
culated light curve is compared with temperatures in- 
ferred from observations in figure 111! We emphasize 
that the calculated blackbody temperatures are recov- 
ered using the luminosity, L(t), and identifying the pho- 
tosphere, iiph(t), with the position where the optical 
depth is t = 2/3, 



■ Calculation 
Observation 



(30) 



Clearly this determination of the blackbody tempera- 
ture is a crude one, and, in fact, the blackbody assump- 
tion mu st generally be con s idered only as an approxi- 
mation (iNakar fe SarillgOlOt iRabinak k Waxmanl l2ml 



IChevalier fc Irwidl201lD~ Therefore, we view the fit be- 
tween calculated and observed temperatures as indicative 
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Fig. 11. — Calculated and observed blackbody temperature of 
SN 2010gx. The calculated temperature (solid black line) was cal- 
culated with the same model as in figure ITDI The observed temper- 
ature (blue circles with error-bars) is taken from Pastorcllo et al. 

poich . 

Considering the simplifications of our approach, we do 
not attempt to find a "best fit" model for the progenitor 
system, but rather demonstrate that a plausible range ex- 
ists. Moreover, we emphasize that the wind parameters 
we choose to fit are not unique. If we fix R w = 2.5 x 10 15 
cm, as implied by the temperature measurement, we have 
a considerable degree of freedom to choose K. The rea- 
son is the weak dependence of the diffusion time on K for 
this R w , as seen in figure [6] The total radiated energy is 
also not strongly effected by an increase of K , as seen in 
figure |H This is because M w « 4irKR w is close to the 
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saturation region of equation (fT0|) and figure [8] More 
specifically, a different model with K = 2 x 10 18 g cm -1 
(double the chosen value), results in a light curve with 
the same peak luminosity as the model in hgurc [10] and 
FWHM wider by only 15%. In addition, we tested a 
model with k = 0.2 cm 2 g _1 , which is appropriate for a 
hydroge n poor CSM, as indic ated by SN 2010gx obser- 
vations (jPastorello et al.ll2010l ). This change has a small 
(less than 10%) effect on the light curve and a negligible 
effect on the temperature, due to the change from inverse 
to linear dependence on k in equation (|17|) (which results 
in a weak dependence at R w ~ Rd)- 

6.2. SN 2006gy 

We take the data for SN 2006gy from ISmith et al.l 
(|2010[) . The measured luminosity and temperature im- 
ply a blackbody radius ~ 4.5 x 10 15 cm at peak luminos- 
ity. If we do not consider the temperature constraint, a 
model which fits the light curve can be found. The model 
with R w = 6.5 x 10 15 cm, K = 0.9 x 10 18 g cm" 1 , and 
E = 5 x 10 51 erg, which reasonably fits the light curve 
except for the tail, is plotted in figure [12] (solid black 
line). We shall distinguish this model as model A. The 
wind mass in this model is M w ~ 35 M^, which is similar 
to the estimates of previous works (IWooslev et al.l [20071 : 
ISmith et al.ll2010t: IChevalier fe Irwin! 120 111) . 

Model A may fit the light curve, but since R w in this 
model is larger than the implied blackbody radius, it de- 
viates from the temperature measurements, as can bee 
seen in figure [T3l We find that a model in a different 
range of the progenitor parameters can be chosen to re- 
produce the observed temperatures (assuming a black- 
body emission) and the peak luminosity, but at the ex- 
pense of generating a light curve which is too narrow. 
One such model, distinguished as model B, is also plotted 
in figures [T2l and [T3l In this model we set R w = 4 x 10 15 
cm, K = 10 18 g cm" 1 , and £ = 3.7x 10 51 erg. 

The light curve is narrow due to the strong dependence 
of the diffusion time on R w , which allows little freedom 
(see § 15. 3|) . As a result, we cannot find a single model 
which fits well both the light curve and the measured 
temperature. This is easily understood by considering 
the analysis in § 15.21 If we use the wind radius fixed by 
the observed peak luminosity and temperature, then for 
a given progenitor star only two free parameters remain: 
K and E. Model B (which assumes the observation- 
ally inferred R w ) is adjusted to reproduce the observed 
peak luminosity, but it results in a narrow light curve, 
yielding a total radiated energy which is too low. The 
total radiated energy cannot be increased by increas- 
ing K, because at K = 10 18 g cm -1 the efficiency of 
converting the explosion energy is already close to the 
asymptotic value (see figure [5] we note that in this figure 
Rw = 2.5 x 10 15 cm; for R. w = 4x 10 15 cm, asymptotic 
efficiency is reached for even smaller values of K, due 
to the larger wind mass). Stipulating a larger explosion 
energy E can account for the total radiated energy, of 
course, but it is not a solution, since it leads to even nar- 
rower light curves, due to larger expansion velocities (see 
figure [7]) . In theory, the desired effect can be obtained 
by increasing both the total energy and K (the latter 
compensating for the narrowing of the light curve, since 
it generates a more massive wind). However, the depen- 
dence of the diffusion time on K is weak for R w in the 



range of interest (see figure [6} , and we determine that 
in order to reproduce the observations, an unrealistically 
heavy wind mass (hundreds of solar masses) is required. 

We note that the fit with observations cannot be im- 
proved by stipulating a larger mass for the progenitor 
star. While an increasing M* leads to a decrease in u s h 
and thus to wider light curves, it works to reduce the ef- 
ficiency of converting the explosion energy to radiation. 
We found that even tripling M* (together with M w , to 
keep the efficiency) , which again results in unrealistically 
large masses, does not widen the light curve enough. 

In essence, it appears difficult to reconcile both the 
observed light curve and temperature of SN 2006gy 
with a steady wind model, since the light curve time 
duration implies a wind radius R w that is differ- 
ent from the one inferred by the photosphere ra- 
dius. T he reason may be the blackbody interpreta- 
tion (see iNakar fc SarilfgOlOl: iRabinak fc Waxmanl 120111: 
IChevalier fc Irwinl l201l[ T which effects the calculated 
temperature, an inaccurate assumption of full ionization 
(which effects the opacity, and therefore the photosphere 
radius), or an indication that a diff erent model, perhaps 
with a different CSM profile (see iMoriva fc Tominagal 
l2012t for example), is required. 

Other, similar, models for SN 2006gy have been sug- 
gested in previous works. IChevalier fc Irwinl (|2011f ) sug- 
gested a steady wind model with estimated parame- 
ters R w ~ 10 16 cm, E ~ 3 x 10 51 erg and K ~ 
0.5 x 10 18 g cm" 1 , resulting in M w ~ 3OM . These pa- 
rameters are similar to the ones we adjusted. Note that 
the large R w in this model resul ts in a low bla c kbody 
temperature, as in our model A. IWooslev et al.l (|2007f ) 
considered a pulsational pair instability scenario, where 
a shell of <~ 30M Q was ejected to a radius of ~ 10 16 
cm prior to the explosion. In this model, the photo- 
sphere radius remains relatively low (this is due to the 
shell density profile and the large radii, which allow for 
a lower optical depth for the same CSM mass) but ap- 
parently n ot low enough to fit to the temperature mea- 
surement. IChatzopoulos et al.1 (|2012h try to explain SN 
2006gy with a semi- analytical model based on the self- 
similar solutions of IChevalie"r] (|1982l ). However, their 
model has a low wind mass (5M Q compared to 40M Q 
ejecta mass), and a small CSM radius (2.5 x 10 15 cm, 
albeit with a different power law), which does not al- 
low to recover the typic al time scale in SN 2006gy. This 
discrepan cy is noted by IMoriva et al.l (|2012l ) . who calcu- 
lated the IChatzopoulos et al.l ( 20121 ) model numerically. 
IMoriva et al.l (]2012t ) suggest their own models for SN 
2006gy, which include interaction with CSM shells which 
are not due to steady winds, but rather a finite shell situ- 
ated at a standoff distance form the star (see their models 
D2 and Fl). The CSM in these models extends to larger 
radii compared to our models (1 x 10 16 — 2 x 10 16 cm), 
and contain less mass (15M© — 18M©), which is com- 
pensated by a greater ejecta energy (10 x 10 51 erg). A 
more detailed comparison with these different models is 
beyond the scope of the current work, but we emphasize 
that our model A fits the light curve with at least the 
same quality as the other numeric a l mod els. One of the 
important aspects of IMoriva et al.1 ()2012l ) is the compar- 
ison of effective and color temperatures (their figure 8). 
The significant difference between the two temperatures 
may be the reason why a simple blackbody approxima- 
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tion cannot account for the observed temperature. 
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in figure 1 141 and the (calculated and observed) temper- 
atures are plotted in figure [T5] The fit of the model to 
the observations is marginal, and suffers from problems 
similar to the fit to SN 2006gy, but more mildly. The 
model is a bit too narrow and the calculated temper- 
ature is a bit too low. This implies some discrepancy 
between the wind radius imposed by the diffusion time 
scale and the blackbody radius imposed by the temper- 
ature and maximum luminosity. We note that that as in 
§ 16.11 the results are not sensitive to the value of k, which 
is s et to be smaller due to lack of hydrogen (as indicated 
by iQuimbv et al.l 120 111) . The discrepancy is again dif- 
ficult to resolve by a different choice of parameters in 
the context of our model, but is more likely to be solved 
if we allow for corrections to the blackbody assumption 
and interpretation of observations. As in the case of SN 
2006gy, if we do not constrain R w by the implied black- 
body radius (because of the non equilibrium conditions 
for example), it is easy to fit the light curve better by 
increasing R w . 



Fig. 12. — Calculated and observed light curves of SN 2006gy. 
Model A (solid black line) was calculated with the parameters 
Rw = 6.5 x 10 15 cm, K = 0.9 x 10 18 g cm" 1 , and E = 5 x 10 51 
erg. Model B (dashed red line) was calculated with the parameters 



Ru 



4 x 10 15 



K 



10 18 g cm" 1 , and E = 3. 7 X 10 51 erg 



The o bserved light curve (blue circles) is taken from Smith ct al. 
(2010), with bolometric correction. 
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Fig. 14. — Calculated and observed light curves of SN 2005ap. 
The calculated light curve (solid black line) was calculated with 
the parameters R w = 2.5 X 10 15 cm, K = 1.5 X 10 18 g cm -1 , and 
E = 2.8 X 10 51 erg. The obs erved light curve (blue circles) is taken 
from lPastorello et alj (f20Toh . with zero bolometric correction. 
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Fig. 13. — Calculated and observed blackbody temperature of 
SN 2006gy. The models (A - solid black line, B - dashed red line) 
were calculated with the same parameters as in figure 1121 The 
observed temperat ure (blue circles with error-bars) is taken from 
ISmith et aLl pTioh 



6.3. SN 2005ap 

The light cur v e dat a for SN 2005ap is taken from 
iPastorello et al.l (12010D an d the temperature is taken 
from IQuimbv et al.l ( 20071 ). The measured luminosity 
and temperature imply a blackbody radius ~ 2.5 x 10 15 
cm at peak luminosity. We choose a model with R w — 
2.5 xlO 15 cm, K = 1.5 xlO 18 g cm" 1 , and E = 2.8 xlO 51 
erg. The calculated and observed light curves are plotted 



7. CONCLUSIONS AND DISCUSSION 

In this work we considered a scenario for superlumi- 
nous supernovae (SLSNe) based on shock breakout from 
the interaction of an ejecta from an exploding star and a 
massive envelope of circumstellar material, presumably 
emitted from the star prior to the explosion. We focused 
on a steady wind model for this CSM, which dictates 
ap~ r~ 2 density profile, and on massive winds, which 
have a total mass M w comparable with the ejected mass, 
Mcj . The latter assumption allows for shock breakout in 
which the entire ejecta participates, leading to an effi- 
cient conversion of the explosion energy to bolometric 
luminosity. This efficiency is required to explain the to- 
tal energy observed in some SLSNe without resorting to 
high energy explosion mechanisms. 

Our approach combines analytical limits and numer- 
ical hydro-diffusion calculations of the bolometric light 
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Fig. 15. — Calculated and observed blackbody temperature of 
SN 2005ap. The calculated temperature ( solid black line) was cal- 
culated with the same model as in figure | 14 | The obser ved tem- 
perature (blue cir cles) is taken fromlQuimbv et al.l 11200 71 ). No es- 
timate is given in Quimbv ct al. (2007]) for the uncertainty in the 
two temperatures in the graph (only a general estimate of 16,000K- 
20.000K). 



curve. Thus we improve upon previous works which con- 
sidered a massive wind in the context of SLSNe with 
order of mag nitud e estimates (ISmith fc McCravl 120071: 
ISmith et all 120101 IChevalier fe Irwinl 120 111) . We note 
that the numerical calculations are absolutely necessary 
for producing reliable results in a M w ps M e j scenario, 
since the analytical limits cannot be applied. On the 
other hand, our numerical model includes only the fun- 
damental physics of the problem, and so the principle 
trends are easily understood. Our calculations are in 
agre ement with more complex numerical models (see, 
e.g., Ilvloriva et a!]|201lD . 

Specifically, we related the main features of the ob- 
served light curve to the parameters of the progenitor 
system - the star and the wind. As expected, we explic- 
itly hnd that a large wind mass does allow to efficiently 
convert the energy of the explosion into a bolometric light 
curve, and that this efficiency naturally saturates when 



M„ 



> M e j. However, we demonstrate that in this sce- 



nario the width of the light curve is strongly dependent 
on the cutoff distance of the wind, R w . Quantitatively, 
we show that in order to recover a time scale of tens of 
days seen in the most energetic SLSNe, the wind outer 
radius must be of order a few 10 15 cm. This result is 
consistent with the scenario of a steady mass loss of a 
few 10~ 2 M f 7)yr~ 1 during the last 10 00 years of the star's 
evolution (jQuataert fc Shiodell2012l ). 

We applied our model to three SLSNe, namely SN 
2010gx, SN 2006gy and SN 2005ap, which exhibited lu- 
minosities as high as several 10 44 ergs with time scales 



of tens of days. We find that the light curves of these 
objects can be understood in terms of a supernova in 
a heavy wind scenario, where the progenitor star has 
a mass of Af* « 15Af Q and the wind mass is in the 
range 15M Q — 35M Q . The efficient conversion of the 
ejecta kinetic energy to radiation in this scenario en- 
ables us to reproduce the light curves with energies in 
the standard range of 2 x 10 51 erg — 5 x 10 51 erg, and, 
as mentioned above, with winds extending to radii of 
2.5 x 10 15 cm — 6.5 x 10 15 cm. Assuming blackbody 
emission, our model reproduces the observations in SN 
2010gx, but underestimates the temperatures measured 
in SN 2006gy, and to some extent in SN 2005ap as well. 
A similar trend was found in the other works mentioned 
above as well, and may be the result of a frequency de- 
pendent (rather than constant) opacity, differences be- 
tween the effective and color temperatures, or both. 

Wc note that a separate subclass of SLSNe, where the 
peak luminosity is high but the time scales are shorter , 
such as SN 2008es (jGezari et al.|[2009l : iMiller et al.ll2009( ) 
may also be explainable with a steady wind model. In 
this case the wind mass must be significantly lower than 
the ejecta mass, so the initial wind optical depth may be 
small. Wc plan to investigate such systems numerically 
in future work. 

We do not explicitly compare our results with an alter- 
native, but similar, scenario, where the ejecta interacts 
with a single shell, situated at some distance from the 
star, presumably a result of a short period of enhanced 
mass loss. Clearly, a massive shell can be equally ef- 
ficient in convertin g the explosion energy to radiation 
(|Moriva et alj [20T2T ) . Differentiating between both sce- 
narios must involve a wider parameter survey (especially 
since in the shell model its position and density struc- 
ture can be assumed independently), and probably some 
additional physics in the numerical model as well. 

Another issue which requires further work is a non- 
thermal component from shock brea kout through the 
thick wind. Recently, several authors dKatz et al.l 120111 : 
IChevalier fc Irwinl 120121 ISvirski et al.ll2012t ) pointed out 
that during breakout the shock is likely to become col- 
lisionless, hence creating a higher energy, non-thermal 
component in the spectrum. These works vary consid- 
erably concerning their estimates regarding the fraction 
of the total non-thermal energy eventually emitted from 
breakout, but obviously, if this fraction is sizable, there 
will be some impact on the bolometric light curve. 
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APPENDIX 

COMPARISON TO PLANAR SHOCK BREAKOUT SOLUTION 



ISakurail (|1960l ) investigated the problem of a shock wave propagating through a non-uniform medium (ideal gas) of 
decreasing density and reaches a boundary where the density vanishes. For a planar initial density profile of p(x) cx x n , 
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with x the distance from the boundary, a self similar solution exists. ISapir et al.l ([201 ll ) studied an extension to the 
problem which includes radiative flux, in the diffusion a pproxima t ion. At large optical depth from the boundary, 
the purely hydrodynamical (without diffusion) solutio n of ISakurail (Il960t ) is applicable, but at r ~ Z?^ 1 , diffusion is 
significant and must be taken into account (see Sj3|). ISapir et al.l (|2011h present a self similar numerical solution to 
the problem, assuming constant opacity and radiation dominated gas. This problem, which can describe the planar 
phase of shock breakout in the absence of wind served as one of several code checks for our program (see § E|) . 

We use planar geometry, insert the appr opriate dens i ty pro file (with n = 3), and keep only the radiation terms of 
the EOS (see § [4]) for the comparison with ISapir et al.l (|2011l) . Our calculation was carried out by moving the inner 
boundary as a piston at a constant velocity. In addition, we deposited thermal energy in the inner most cell as an 
initial condition. The shock wave which arises in these conditions converges to the self similar solution (|Sakurailll960l) 
as it propagates through the gas. In figure IA16I we dem onstrate the exact fit that we find between our calculated light 
curve and the self similar solution of ISapir et al.1 ( 20111 ). Our light curve was normalized to breakout point (r = f3~^) 
related values to and Cq — pqVq (for details see ISapir et al.ll201ll) . 



10 u 



10 



10 



10 



Our Code 

O Sapir et al. 




20 40 60 

(* ~ *peak)/*o 



100 




5 

(* - *peak)/*o 



10 



Fig. A16. — Normalized (sec Sapir ct al. 2011) emitted energy flux as a function of normalized time relative to time of peak emitted 
energy flu x, for a den s ity pr ofile with n = 3. The solid black line is the curve calculated by our code, and the blue circles are taken from 
table 3 of ISapir et all i201lD . The right figure is a zoom of the left figure in a shorter time span near the peak. 
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